Exact Mapping of the 2+1 Dirac Oscillator onto the Jaynes-Cummings Model: 

Ion- Trap Experimental Proposal 



o 
o 

(N 

Oh- 
< 
OO 



■ 

9 L,: 



> 
in 

m 

(N 
O 

o 



x 



A. Bermudcz 1 , M. A. Martin-Dclgado 1 and E. Solano 2 ' 3 
1 Departamento de Fisica Teorica I, Universidad Complutense, 28040 Madrid, Spain 
2 Physics Department, ASC, and CeNS, Ludwig-Maximilians-Universitdt, Theresienstrasse 37, 80333 Munich, Germany 
3 Seccion Fisica, Departamento de Ciencias, Pontificia Universidad Catolica del Peru, Apartado Postal 1761, Lima, Peru 

We study the dynamics of the 2+1 Dirac oscillator exactly and find spin oscillations due to a 
Zitterbewegung of purely relativistic origin. We find an exact mapping of this quantum-relativistic 
system onto a Jaynes-Cummings model, describing the interaction of a two-level atom with a quan- 
tized single-mode field. This equivalence allows us to map a series of quantum optical phenomena 
onto the relativistic oscillator, and viceversa. We make a realistic experimental proposal, at reach 
with current technology, for studying the equivalence of both models using a single trapped ion. 
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Current technology has allowed the implementation of 
the paradigmatic nonrelativistic quantum harmonic os- 
cillator in a single trapped ion [l[ , one of the most funda- 
mental toy models in any quantum mechanical textbook. 
However, its relativistic version, the so-called Dirac oscil- 
lator 0, [H , remains still far from any possible experimen- 
tal consideration for different fundamental and technical 
reasons. We will show here that available experimen- 
tal tools may allow the implementation of the relativistic 
Dirac oscillator in a single nonrelativistic trapped ion. 

The Dirac oscillator was introduced as an instance of 
a relativistic wave equation such that its nonrelativistic 
limit leads to the well-known Schrddinger equation for 
the harmonic oscillator. This is achieved by introducing 
the following coupling in the Dirac equation 



cctj (p* — imfiur 3 ^ + fimc 2 
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(i) 



where |^/) is the Dirac 4-component bispinor correspond- 
ing to a quantum relativistic spin-i particle, like the elec- 
tron, c is the speed of light, m is the particle rest mass, 
and a.j, /3, arc the Dirac matrices in the standard repre- 
sentation. The interacting Hamiltonian is linear in both 
momentum p } and position r 3 , j = x,y,z, and lo turns 
out to be the harmonic oscillator frequency. Remark that 
when uj = we recover the standard Dirac equation [ij]. 
The Dirac oscillator looks like a particular gauge trans- 
formation p — > p |A that is linear in position, but the 
presence of the i and the f3 matrix makes a crucial differ- 
ence. Demanding the correct energy-momentum relation 
for a relativistic free particle E = y / p 2 c 2 + m 2 c 4 , these 
matrices are 4x4 dimensional and must obey a Clifford 
algebra given by the anticommutation relations 



ajCXk + akCtj 
otj/3 + /3a j 
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There has been a growing interest in simulating quan- 
tum relativistic effects in other physical systems, such as 



black hole evaporation in Bose-Einstein condensates Q 
and the Unruh effect in an ion chain Another as- 
tonishing relativistic prediction is the Zitterbewegung [H , 
a helicoidal motion realized by the average position of a 
relativistic fermion, which has been discussed in the con- 
text of condensed matter systems Q and the free-particle 
Dirac equation in a single ion 0. 

Here, we shall be concerned with the Dirac oscillator 
model in 2+1 dimensions, since it is in this setting where 
we can establish a precise equivalence with the Jaynes- 
Cummings (JC) model Q. In two spatial dimensions, 
the solution to the Clifford algebra ([2]) is given by the 
2x2 Pauli matrices: a x = a x ,a y = a y ,(3 = a z . In 
this case, Y$>) can be described by a 2-component spinor 
which mixes spin up and down components with positive 
and negative energies. In particular, the Dirac oscillator 
model now takes the form 
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In this paper, we shall provide the complete (eigenstates 
and energies) and exact solution of the 2D Dirac oscilla- 
tor in order to study its relativistic dynamics, where cer- 
tain collapses and revivals in the spin degree of freedom 
appear as a consequence of Zitterbewegung. In addition, 
we derive an exact mapping of the 2+1 Dirac oscilla- 
tor onto the JC model, an archetypical quantum optical 
system. Furthermore, we propose the simulation of this 
relativistic dynamics in a single trapped ion, a physical 
setup possessing outstanding coherence features. 



Considering the spinor 



tion ([3|) becomes a set of coupled equations 
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In order to find the solutions, it is convenient to introduce 
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the following chiral creation and annihilation operators 



a r 
ft, 



-±(a x -ia y ), 4 : = 7f(4 + iat), 
^(a^+iay), oj := ^=(4 - iaj), 



V2 V 



where a x , aj., a y , a y , are the usual annihilation and 



creation operators of the harmonic oscillator a 
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and A = yj h/muj represents the 
ground state oscillator width. The orbital angular mo- 
mentum may also be expressed as 



L z = h(a].a r — a\ai), 



(6) 



which leads to a physical interpretation of a\ and a\. 
These operators create a right or left quantum of angular 
momentum, respectively, and are known hence as circu- 
lar creation-annihilation operators. Equations (|4|) can be 
rewritten in the language of these circular operators 



l^i 
life 
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a l I'M, 

ai \ipi), 



(7) 



where £ := hu/mc 2 controls the nonrelativistic limit. In 
order to find the energy spectrum we shall solve the asso- 
ciated Klein-Gordon equation, which can be derived from 
Eqs. ([7]) as follows 



(E 2 - m 2 c 4 )|^) = 4m 2 c 4 £ a\ai |V>i>, 
(E 2 - m 2 c 4 )|^ 2 ) = 4m 2 c 4 £ (1 + oJo,)|^>- 
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These equations can be simultaneously diagonalizcd writ- 
ing the spinor in terms of the left chiral quanta basis 



f 



vac) 



where m = 0, 1, ... The energies can be expressed as 

(^ 2 i -m 2 c 4 )|n;)=4m 2 cVi |m>, 
(S 2 ; -m 2 c 4 )K)=4m 2 c 4 e (l+n{)|n{). 



(9) 
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Since both components l^i) and ^2) belong to the same 
solution, the energies must be the same E n i = E ni . This 
physical requirement sets up a constraint on the quantum 
numbers n; =: n[ + 1. Note that, following ©, the state 
\m) corresponds to a negative angular momentum. The 
energy spectrum can be described as follows 



E = ±E ni = ±mc 2 y / l + 4:£ni. 



(11) 



To find the corresponding cigenstates, we go back to 
Eq. (|7|) , and after normalization we arrive at the expres- 
sion for the positive and negative energy eigenstates 



±E ni ) = 
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where the quantum number is now restricted to ni = 
1,2, ... In this way, we have solved the two-dimensional 
Dirac oscillator describing the energy spectrum and the 
eigenstates in terms of circular quanta. The distinction 
between Dirac and Klein-Gordon eigenstates is an im- 
portant point in order to understand the dynamics of the 
2+1 Dirac oscillator and its realization in an ion trap. 

The eigenstates of the 2D Dirac oscillator can be 
expressed more transparently in terms of 2-component 
Pauli spinors |%t) and 1X4 ) 



E ni ) 



*niK)IXT> -i/?n>i - l)|Xi), 
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(13) 



are real. 



From this expression we observe that the energy eigen- 
states present entanglement between the orbital and spin 
degrees of freedom. This property is extremely important 
since the following initial state 

|*(0)) := \ ni -l)\ Xl )=iPm\+E ni )-ia ni \-E ni ) (14) 

superposes states with positive and negative energies, 
and this is the fundamental ingredient that leads to Zit- 
terbewegung in relativistic quantum dynamics. This phe- 
nomenon, due to the interference of positive and negative 
energies, has never been observed experimentally. The 
reason is that the amplitude of these rapid oscillations 
lies below the Compton wavelength, where pair creation 
is allowed, and the one-particle interpretation falls down. 

Now, the evolution of this initial state can be expressed 
in the energy basis as 



|*(t)>=i/3 n ,|+^ l )e- u, **«* 



E n ,)e 



where 



E n , mc 2 1 — — 

~f = — 



(15) 



(16) 



describes the frequency of oscillations. Writing this 
evolved state in the language of Pauli spinors, 




sinw„,t \m - 1)1X4.) 



smu ni t \ni)\xi), 



(17) 



we observe an oscillatory dynamics between |ti;)|xt) an( l 
\ni — 1)1X4.)- I"! 10 hritial state, \ni — l)|xj,), which has 
spin-down and n\ — 1 quanta of left orbital angular mo- 
mentum, evolves exchanging a quantum of angular mo- 
mentum from the spin to the orbital motion. 

The dynamics described in (|17|) is completely similar 
to the atomic Rabi oscillations occurring in the Jaynes- 
Cummings model, though arising from a completely dif- 
ferent reason. Whereas the Rabi oscillations in the 
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Jaynes-Cummings model are caused by the interaction of 
a quantized electromagnetic field with a two-level atom, 
the relativistic oscillations are caused by the interference 
of positive and negative energy states and therefore con- 
stitute a clear signature of Zitterbewegung 0] • 

To clarify this issue further, we calculate the time evo- 
lution of the following physical observables, that catch 
the full essence of the system dynamics, 



(L z ) t = -(nj 
(S*)t = 



i)h 



^^ftsin 2 u ni t, 



(18) 



where J z = L z + S z stands for the ^-component of the 
total angular momentum. The latter relations describe 
a certain oscillation in the spin and orbital angular mo- 
mentum, while the total angular momentum is conserved 
due to the existent invariance under rotations around the 
z-axis. It is important to highlight that these oscillations 
have a pure relativistic nature. In the nonrelativistic 
limit ^ < 1, these oscillations become vanishingly small 

{L z ) t = -(nj - l)h-^mhsm 2 a ni t + 0(e), hq) 
(S z ) t = -| +4£n l h S m 2 n ni t + 0(e), 1 ' 

where Q, ni := mc 2 (l + 2(ji{)/h stands for the oscillation 
frequency in the nonrelativistic limit. In this limit the 
negative energy components are negligible and therefore 
the Zitterbewung disappears. 

The results discussed so far allow a precise mapping 
between two seemingly unrelated models: the Jaynes- 
Cummings model of Quantum Optics and the 2D Dirac 
oscillator. Starting from Eq. we may write the Dirac 
oscillator Hamiltonian as 



H = 2imc 2 ^l (a] \ip 2 )(A I - a/IV'iXV'af) + mcV 



h(g<7 a\ + g*a + a{) + mc 2 a z , 



(20) 



where er + , a~ , are the spin raising and lowering opera- 
tors, and g := 2\mc 2 ^f£ l /h is the coupling strength be- 
tween orbital and spin degrees of freedom. In Quantum 
Optics, this Hamiltonian describes a Jaynes-Cummings 
interaction, that has been studied in cavity QED and 
trapped ions [l], [l(|, among others. Within this novel 
perspective, the electron spin can be associated with a 
two-level atom, and the orbital circular quanta with the 
ion quanta of vibration, i.e., phonons. As we will see be- 
low, the central result of Eq. ([20]) allows both physical 
systems, the JC model and the 2D Dirac oscillator, to 
exchange a wide range of important applications. 

We will show now how to implement the dynamics of 
Eq. ((3|) in a single ion inside a Paul trap, which was shown 
to follow the dynamics of Eq. (|20|) . The Dirac spinor will 
be described by two metastable internal states, \g) and 
\e), as follows 



while the circular angular momentum modes will be rep- 
resented by two ionic vibrational modes, a x and a y . Cur- 
rent technology allows an overwhelming coherent con- 
trol of ionic internal and external degrees of freedom [lj . 
There, three paradigmatic interactions, the carrier, red- 
, and blue-sideband excitations, can be implemented at 
will, independently or simultaneously [111 ]. For example, 
using appropriately tuned lasers, it is possible to produce 
the following interactions 



H AJC 



a + aU iv + a-"a<e" 



(22) 



where {a,i,a\}, with i = x,y, arc the phonon annihila- 
tion and creation operators in directions x and y, vi are 
the natural trap frequencies, r\i := ki^Jh/2Mvi are the 
associated Lamb-Dicke parameters depending on the ion 
mass M and the wave vector k, 5i and Qi are the ex- 
citation coupling strengths and 0, <p, the red and blue 
sideband phases. Remark that the term h5i<j z , in H? c of 
Eq. (|2"2"1) , stems from a detuned JC excitation. 

A suitable combination of the above introduced exci- 
tations (|22|). with proper couplings and relative phases, 
can reproduce the following Hamiltonian 

H = c [af Vx + afpy] + mujc [afy - afx] + mc 2 af 

(23) 

with of := \g) (e| + \g)(e\, of := -i(\e){g\ - \e){g\), 
<jf e := |e)(e| — |.g)(<?|, and the following parameter corre- 
spondence 



c = a/^OA, 
mc 2 = h5, 
mujc = 



(24) 



where A := Aj is the width of the motional ground state, 
Q := Qi,r/ := 7?;,Vi = x, y. The remarkable equivalence 
of the Dirac oscillator Hamiltonian ((3|) and the interac- 
tion in Eq. (|23|) shows that it is possible to reproduce 
the 2D Dirac oscillator, with all its quantum relativis- 
tic effects, in a controllable quantum system as a single 
trapped ion. 

For the sake of illustration, note that the effective 
terms appearing in Eq. (|23p can be achieved by suitable 
linear combinations of Hf G and ff l AJC in (f2"2"| . 



i = x,S x = 5,4>= 2z,ip= f 
i = y, S y = 0, <f) = , ip = 7r 
i = x,S x = 0,<f> = §, ip = f 
i = y, 5 y = 0, <j) = , ip = 



Px 



V2fir]f2Ao- 
V^hrjnAafpy, 
V2Hri£lA- 1 Oy' e . 
s/2TvqQ.h- 1 o s x 



y- 



W)\e) + \<h)\9) 



(21) 



(25) 

Note that in the trapped ion picture, the important pa- 
rameter £ = 2(r)(l/S) 2 can take on all positive values, 
assuming available experimental parameters: 77 ~ 0.1, 
Q, ~ - 10 6 Hz, and S ~ - 10 6 Hz fl. The ability to ex- 
perimentally tunc these parameters will allow the exper- 
imenter to study otherwise inaccessible physical regimes 
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that entail relativistic and nonrclativistic phenomena. 
For example, the Zitterbewegung is encoded in the spin 
degree of freedom, and we can associate Rabi oscillations 
to the interference of positive and negative energy so- 
lutions. Setting the initial state |0)|xj.) *-» |0)|<?}, the 
internal degree of freedom evolves according to Eq. (fig)) 



h 

2 



ft sin u>it, 



l + 4£ 



(26) 



where u)\ = S\/l + 4£, see Eq. (jl6[) . stands for the fre- 
quency of the Zitterbewegung oscillations and can take 
on a wide variety of measurable values. 

In order to simulate this dynamics in an ion-trap table- 
top experiment, the ion must be cooled down to its vi- 
brational ground state |0), with a current efficiency above 
99% Q. To estimate the observable ([26]) . one can make 
use of the powerful tool called electron shelving, where 



(S z 



[2P e {t) - 1] 



(27) 



can be obtained through the measurement of the proba- 
bility of obtaining the ionic excited state P e (t) with ex- 
traordinary precision. 

Another fundamental result of the JC model which can 
be mapped straightforward to the Dirac oscillator is the 
existence of collapses and revivals in the atomic popu- 
lation, which is claimed to be a direct evidence of the 
quantization of the electromagnetic field. To produce 
this effect an initial state \z)\g) is required, where \z) is 
an initial circular coherent state, 



measured via a similar electron-shelving technique, while 
the observable of Eq. (|30|) can be measured via the map- 
ping of the collective motional state onto the internal 
degree of freedom 

It is worth mentioning that the chiral partner of the 2D 
Dirac oscillator Hamiltonian © can be obtained through 
the substitution lo — > — uj, and consists on right-handed 
quanta. This Hamiltonian presents similar features as 
those discussed above, and can be exactly mapped onto 
an anti-Jaynes-Cummings interaction 



H = h(ga r a + g*a\a + ) + mc 2 cr z , 



(31) 



with similar parameters. It is precisely this chirality 
which allows an exact mapping between the JC, AJC, 
and the lefthanded and righthanded 2D Dirac oscillator. 
This essential property, missing in the 3D case, forbids 
an exact mapping of Eq. ((T|) onto a JC-like Hamiltonian. 

In conclusion, we have demonstrated the exact map- 
ping of the 2+1 Dirac oscillator onto a Jaynes-Cummings 
model, allowing an interplay between relativistic quan- 
tum mechanics and quantum optics. We gave two rel- 
evant examples: the Zitterbewegung and collapse-revival 
dynamics. In addition, we showed that the implementa- 
tion of a 2D Dirac oscillator in a single trapped ion, with 
all analogies and measured observables, is at reach with 
current technology. 
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I*(0)> 



"1=0 



\ni)\g), 



(28) 



with z G C. After an interaction time t, 

|2«!p-M 2 



h f, 4Z(m + l)\z\-"'e 
2 _n [l + 4£(n, + l)]n,! 



sin 2 (cj„ !+ ii). 

(29) 

This expression can be understood as an interference ef- 
fect of terms with different frequencies Lu nt +i leading to 
collapses and revivals. A novel feature of the Dirac oscil- 
lator is the appearance of these collapses and revivals in 
the orbital circular motion of the particle, reflected in 



(L z ) t = -Hz 2 -h , ' A J ' m , sin 2 K,+i* ■ 

(30) 

The generation of an initial circular coherent state will 
require two sequential applications of the technique dc- 
scribed in Ref. [1[ on an initial motional ground state. 
These two operations should be applied with a relative 
phase such that Di(z) = D x (z)D y (—iz) i where Dj(z) = 



j = x,y. The observable of Eq. (|29|) can be 
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